Phase Transitions in Non-extensive Spin Systems 
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The spherical spin model with infinite-range ferromagnetic interactions is investigated analytically 
in the framework of non-extensive thermostatics generalizing the Boltzmann-Gibbs statistical me- 
chanics. We show that for repulsive correlations, a new weak-ferromagnetic phase develops. There is 
a tricritical point separating para, weak-ferro and ferro regimes. The transition from paramagnetic 
to weak-ferromagnetic phase is an unusual first order phase transition in which a discontinuity of 
the averaged order parameter appears, even for finite number of spins. This result puts in a new 
way the question of the stability of critical phenomena with respect to the long-ranged correlations. 
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Discussion of nonextensivity in thermodynamics go 
back to 70 's p0[ . Recently, non-extensive theories have 
become extremely important in several areas of Physics 
Pj . However, there are still very important questions that 
should be addressed. The influence of the non-extensivity 
on phase transitions and their trace in finite systems is 
certainly a field which should be further explored S . In 
the collisions of atomic nuclei or charged atomic clus- 
ters, a highly excited non-equilibrium transient system is 
formed which at a later stage of the reaction approaches 
an equilibrium and splits into many fragments |^,|| . Ob- 
served signatures of criticality in these processes depend 
strongly not only on strong repulsive Coulomb interac- 
tions but also on on the way the excited transient system 
is formed in collisions || . The fragility of those signatures 
manifests the inherently non-extensive character of these 
'critical phenomena' [M. The general problem of the re- 
lation between critical behavior in small non-extensive 
statistical systems and the phase transitions in the ther- 
modynamical limit is investigated here using the Berlin- 
Kac Model (BKM) g] in the framework of the Tsallis 
generalized statistical mechanics (TGSM) Q . 

The TGSM was inspired by earlier works on non- 
extensive thermodynamics, going back to the 60' [flo|| , 
it is based on an alternative definition for the equilib- 
rium entropy of a system whose ith microscopic state 
has probability pi B : 
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and q (entropic index which is determined by the mi- 
croscopic dynamics) defines a particular statistics. In 
the limit q = 1, one obtains the usual Boltzmann-Gibbs 
(BG) formulation of the statistical mechanics. The main 
difference between the BG formulation and the TGSM 
lies in the nonadditivity of the entropy. For two indepen- 
dent subsystems A, B, such that the probability of A + B 
is factorized into : pa+b = PaPb-, the global entropy ver- 
ifies : S q (A + B) =S q (A) + S q (B) + (l-q)S q (A)S q (B) . 



The entropy S q has definite concavity for all values of 
q. In particular, it is always concave for q > 0, and this 
is the case discussed in the present paper. It has been 
shown |l2] that the TGSM retains much of the formal 
structure of the standard theory. Many important prop- 
erties , such as the Legendre transformation structure of 
thermodynamics and the H-theorem (macroscopic time 
irreversibility) have been shown to be ^-invariant. Con- 
sidering the fact that the essence of the Second Law of 
Thermodynamics is concavity (see Ref. 14 1), the men- 
tioned properties of this entropy allow us to say that 
there are no problems with this law in TGSM. 

The TGSM is relevant if the effective microscopic in- 
teractions are long-ranged and/or the effective micro- 
scopic memory is long-ranged and/or the geometry of the 
system is fractal [p| JT^Jl5[ |. In the superadditive regime 
(1 — q > 0), independent subsystems A and B will tend 
to join together increasing in this way the entropy of 
the whole system. On the contrary, in the subadditive 
regime (1 — q < 0), the system increases its entropy by 
fragmenting into separate subsystems. These ideas are 
in agreement with the results of Landsberg et al. [jll 16 
The subadditivity of entropy is expected when long-range 
repulsive interactions and correlations and/or long-range 
memory effects are present in the system |l3|]l6| 1. This 
is in particular the relevant limit of the fragmentation 
of electrically charged (off-equilibrium) systems, such as 
formed in the collisions of atomic nuclei or sodium clus- 
ters. 

The BKM was introduced as an approximation of Ising 
model. In the BKM, the individual spins are taken as 
continuous 3-dimensional variables, but the rigid con- 
straint on each local spin variable Si (where Si = ±1/2) 
in Ising model, is relaxed and replaced by an overall 
spherical constraint : 
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(2) 



for N spins. In the present work, we consider the infmitc- 



1 



range version of the BKM where all pairs of spins interact 
with the same strength, according to the Hamiltonian : 



H — \ Si Sj 

TV £—1 8 



(3) 



The sum in (^) runs over all different pairs of spin indices 
(i,j). The constant term — J/8 is added to fix the ori- 
gin of the energies at the vanishing total magnetization 
• Si & = 0. The normalization of the ferromagnetic 
(J > 0) exchange energy as J/N, ensures finiteness of 
the energy per spin in an infinite system. The value of 
this energy should depend on the volume of the system. 

One defines the magnetization per site by : r\ = 
N~ x ^2 li S{ . Because of the inequality : {J2Si) 2 < 
Nj^Sf, the value of r\ is contained in between —1/2 
and +1/2. The phase of the system will be characterized 



by the average value of r\ 



2 . 



-< T) > , where the 



brackets < • • • > denote the ensemble average at a given 
temperature. The square root of m 2 plays the role of the 
order parameter. 

The Hamiltonian of BKM can be expressed as a func- 
tion of rj : 
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(4) 



and this allows us to obtain the degeneracy factor of the 
macroscopic states defined by the fixed number of spins 
N and the fixed energy E. The number of such states 
in the ./V-dimensional spin-space is the volume of the in- 
tersection of the hypersphere (|^) and of the hyperplane : 
Y.iSi = {-2NE/J) 1 / 2 . This intersection is a sphere of 
radius (N/4+2E/ J) 1 / 2 in the (N—l)— dimensional space 
and its volume is : g N (E) = a N (l/4 + 2E/JN) N / 2 ~ 1 , 
where a at is a numerical factor which depends only on 
the number of spins (system mass). Knowledge of this 
statistical weight, in addition to the energy given by (Q) 
allows us to compute all thermodynamic quantities of the 
equilibrated system. 

In BG statistical mechanics, the partition function is 
written in terms of the redistribution as : 

,1/2 

Z = A rj exp(-Nf(f]))dn (5) 
Jo 

The constant A is : 2 5 ~ N JNa N exp(-f3 J/8) and : 

f( V ) = -PJ V 2 - (1 - 2/A0M1 - 4ry 2 )/2 

is the free energy per spin. In the limit N — > oo, this 
function does not depend on N which ensures the exis- 
tence of the thermodynamic limit. Since this free energy 
is analytical in the variable rj : f(rj) — [4(1 — 2/N) — 

l3J]rf/2 + 8(1 - 2/N)?f H , with a positive ry 4 -term 

and a change in sign of the ry 2 -term at the pseudo-critical 
point (j3J) c> N = 4(1 — 2/N), the infinite system under- 
goes a second order phase transition at the critical tem- 
perature (/3J)c = 4. The critical exponents are given by 



the regular Landau-Ginzburg mean-field theory for the 
magnetic systems. 

In the TGSM, the equilibrium energy distribution is 
given by M : 



PN,q{E) 



(l-(3(l-q)E) 



1/(1-9) 



(6) 



if 1 — /3(1 — q)E > 0, and otherwise. The partition 
function Z q is defined by the normalization oipN,q ■ 



J g N {E) PN , q {E)dE = 1 



(7) 



All ensemble averages are then defined with respect to the 
properly normalized probability : p q N q (E)g^(E) |fl7| . 
For example, the averaged order parameter is given by : 



Jr, 2 (E)p%^E)g N (E)dE 
m* =< i] > q =- - 



J p%JE)g N (E)dE 



(8) 



Replacing the integration variable E by rj, the value of 
Tfiq is : 



m 2 q = h/h , 
with the integrals / s defined by : 

-i -N/k-1 



(9) 
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lo 
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and x = min(l/2, ^/2//3Jk) because of (||) and \r)\ < 
1/2. The fundamental parameter of the non-extensive 
generalization of the BKM is : 



K = N(q-1) 



(11) 



instead of q — 1. This parameter, which will be called 
hereafter the Out-of-Extensivity (OE) parameter, allows 
us to define the thermodynamic limit of the system in a 
natural way. Let us compare systems of different masses 
N such that the OE parameter ( |TT| ) remains constant. In 
the relations (^|) and (|I(i|), one can recognize the equiva- 
lent of a free energy per site : 



f q (rj) = ln(l - f3J KV 2 /2)/ K - ln(l - 47? 2 )/2 



(12) 



which is independent of N when k is a constant. Exis- 
tence of this free energy per spin implies in turn existence 
of the thermodynamic limit in the system. For a fixed 
value of the OE parameter, the interesting values of q 
approach 1 as the number of spins increases. 

When the number of spins N is finite, the integrals I s 
can be expressed by the hypergeometric functions. At 
high temperatures : (3 J < 8/k, the upper bound for 
rj in ( p~0[ ) is 1/2. The integrand in this case is always 
finite, and m q is a positive decreasing function of the 
temperature. Moreover, if : k > 2/(1 — 2/N), then for 



2 



all these temperatures one has : (3 J < 4(1 — 2/N) and 
m q oc 1/v^V, similarly to the standard paramagnetic 
phase. In particular, when (3 J — > 8/k by lower values 
then the limiting value is m q — N^ 1 ^ 2 [k/(2(k — 2))] 1 / 2 . 
On the other hand if k < 2/(1 - 2/N), then this 1/ViV- 
behavior of m q holds only when (3J < 4(1 — 2/N). For 
higher values of /3J, the leading behaviour of the or- 
der parameter takes finite values independent of N, and 
tends to 1/2 when [3 J — > 8/k by the lower values. 
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FIG. 1. Plots of the averaged magnetization of the infi- 
nite-range BKM in the TGSM for N = 50 interacting spins 
and k = 2k c , k c , k c /2. k c is the value of the OE parameter for 
which the pseudo-critical point of the standard second order 
phase transition (here ((3J) c ,n = 3.84) and the new critical 
point for the non-extensive BKM ((/3J) c ,q = 8/k) occur at 
the same temperature. The BG limit (k = 0) is shown by the 
dashed line. 

At low temperatures : (3 J > 8/k, the quantities 
[1 — {I3Jk/2)t]' 2 ]~ n / k ~ 1 diverge for finite values of r\ 
smaller than 1/2, and the value of m q is dominated by 
the behavior of the integrand near this divergency. In 
this case : m q = y/2/ (J3Jk) , since the order-parameter 
probability distribution collapses into the two Dirac dis- 
tributions. Moreover, if (3 J tends to 8/k by larger val- 
ues, the limit of m^-values is 1/2. This implies that if 
k > 2/(1 — 2/N), then the averaged order parameter 
m q has a discontinuous jump at the critical temperature 
{(3J) c ,q = 8/k between the small value of order l/y/~N 
and 1/2. Otherwise, if k < 2/(1 — 2/N), then m q is 
continuous but its first temperature derivative becomes 
discontinuous at this critical temperature ((3J) c ,q — 8/k, 
leading to the second order critical phenomenon. Both 
behaviors are exemplified in Fig. 1. The tricritical point 
{((3J)tc = 4, ktc = 2) separates the phase bound- 
ary line in transitions of different nature (discontinu- 
ous/continuous) . One of the continuous phase transitions 
is the simple continuation of the well known paramag- 
netic ferromagnetic transition for the BG statistics, 
while the other one is a new second order phase tran- 
sition (ferromagnetic ^ weak-ferromagnetic) due to the 



presence of the long-range repulsive correlations/memory 
effects. 
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FIG. 2. Phase diagram of the infinite-ranged BKM 
ferromagnetic spin model in the TGSM at the thermo- 
dynamic limit. The bold-face solid curve shows the 
line where the first order phase transition between para- 
magnetic and weak-ferromagnetic phase happens. For 
< k < 2, the system undergoes two distinct second or- 
der phase transitions : paramagnetic<^>ferromagnetic and 
ferromagnetic<=>weak-ferromagnetic. The tricritical point is 
located in (/3J)tc = 4, ktc = 2. 

Since the above limiting temperatures and OE param- 
eters do not depend on the mass N (the number of 
spins), the pattern of different behaviors of the order- 
parameter remains unchanged in the thermodynamic 
limit (N — ► oo) and one may draw the phase diagram 
(3 J — k (k > 0) (see Fig. 2). The first order phase transi- 
tion in the generalized BKM (the bold-face solid line in 
Fig. 2) remains unaltered in small systems, unlike the 
case of regular collective critical phenomena. The ad- 
ditional phase appearing on this diagram is due to the 
long-range non-extensive correlations which tend to dis- 
organize spin coherence. This phase is called the weak- 
ferromagnetic state, because the order parameter in this 
state has a positive value, but this value is lowered by 
the disruptive long-range correlations (see Fig. 1 in the 
case k = k c /2). This is the reason why the magnetiza- 
tion in the weak-ferromagnetic phase decreases when the 
temperature is lowered. 

As we have seen above, the energy at constant volume 
in the BKM can be written as a function of the squared 
order parameter. This allows us to write the averaged 
energy as : U q = —JNm q /2 . The specific heat at con- 
stant mass and constant volume is the derivative of U q 
with respect to the temperature, keeping the mass N and 
the coupling J fixed : 



C a 



{(3Jf N dm\ 



d{fij) 



(13) 



,/From the above discussion, it is clear that there is a 
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finite discontinuity of the specific heat when re < 2/(1 — 
2/N). But, there exists also a true divergence of C q for 



the temperature (f3J) c ,q = 
even if the system is finite. 
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FIG. 3. Plots of the specific heat in the infinite-range 
BKM (TV = 50) for : re = re c (the divergence at 
(/3J)c, 9 = 8/re), and re = re c /2 (the jump at (flJ) c , q = 8/re). 
The BG limit (re = 0) is shown by the dashed line. 

There is some evidence |Q supporting the fact that 
long-range attractive interactions lead to q < 1. We be- 
lieve it is very important to stress that long-range repul- 
sive interactions can lead to non-extensivity with q > 1. 
The fact that the present system undergoes two different 
second-order phase transitions in the region q > 1 is very 
remarkable. On the other hand, in this system the fun- 
damental parameter is N(q — 1) instead of q. This is a 
new result in non-extensive physics. 

In conclusion, we have shown that the spherical spin 
model with long-range correlations/memory effects simu- 
lated in the framework of non-extensive thermostatistics, 
develops a new weak-ferromagnetic phase in a subaddi- 
tive entropy regime. An unusual first order phase transi- 
tion, which exhibits discontinuity of the order parameter 
even in finite systems, separates this phase from the stan- 
dard paramagnetic phase. On the phase boundary line 
in the plane (3 J — re, we have found the tricritical point 
separating the nature (discontinuous/continuous) of the 
transition. Above a critical value of the OE-parameter, 
the spin system freezes into the weakly ordered (weak- 
ferromagnetic) phase, passing through the phase bound- 
ary in the discontinuous transition. For a fixed value of 
the OE-parameter, the thermodynamic limit (N — > oo, 
q — 1 — > 0+) of the spin system differs from its BG limit. 
The non-extensivity shields the system from the continu- 
ous disorder <J=^order phase transition and suggests that, 
in general, the second order critical phenomena may be 
unstable in the presence of long-range repulsive corre- 
lations. This result puts in a different perspective the 
discussion of a 'critical behavior' in collisions of atomic 
nuclei or atomic clusters, showing that the observed sig- 



nals correspond possibly to a different limiting behavior 
than previously supposed. 

Acknowledgements 

We thank K.K. Gudima for stimulating discussions. 



[1] A. Renyi, Probability theory (North-Holland, Amster- 
dam, 1970); Z. Daroczy, Information and Control 16, 
36 (1970); L. Galgani and A. Scotti, Pure and Ap- 
plied Chemistry 22, 229 (1970); P. T. Landsberg, Pure 
and Applied Chemistry 22, 215 (1970); J. Aczel and Z. 
Daroczy, On measures of information and their charac- 
terization (Academic Press, New York, 1975); 

[2] R. Salazar and R.Toral, Phys. Rev. Lett. 83, 4233 (1999); 

F. Tamarit and C. Anteneodo, Phys. Rev. Lett. 84, 208 
(2000); M Buiatti, P. Grigolini, and A. Montagnini, Phys. 
Rev. Lett. 82, 3383 (1999); D. B. Walton and J. Rafelski, 
Phys. Rev. Lett. 84, 31 (2000); 

[3] F. D. Nobre and C. Tsallis, Physica A 213, 337 (1995) 

[4] J. P. Bondorf et al, Phys. Rep. 257, 133 (1995). 

[5] F. Chandezon et al, Phys. Rev. Lett. 74, 3784 (1995); 

G. Seifert, R. Gutierrez and R. Schmidt, Phys. Lett. A 
211, 357 (1996). 

[6] J. Pochodzalla et al, Phys. Rev. C 55, 2991 (1997); Y.G. 
Ma et al, Phys. Lett. B 390, 41 (1997); Z. Majka et al, 
Phys. Rev. C 55, 2991 (1997); M.D'Agostino et al, Nucl. 
Phys. A 650, 329 (1999). 
[7] K.K. Gudima et al, Phys. Rev. Lett. 85, 4691 (2000). 
[8] T.H. Berlin and M. Kac, Phys. Rev. 86, 821 (1952). 
[9] C. Tsallis, J. Stat. Phys. 52, 479 (1988). 

[10] A. Renyi, Probability theory (North-Holland, Amster- 
dam, 1970); Z. Daroczy, Information and Control 16, 
36 (1970); L. Galgani and A. Scotti, Pure and Ap- 
plied Chemistry 22, 229 (1970); P. T. Landsberg, Pure 
and Applied Chemistry 22, 215 (1970); J. Aczel and Z. 
Daroczy, On measures of information and their charac- 
terization (Academic Press, New York, 1975); 

[11] P. T. Landsberg and D. Tranah, Collective Phenomena 
3, 73 (1980); D. Tranah and P. T. Landsberg, Collective 
Phenomena 3, 81 (1980). 

[12] A. M. Mariz, Phys. Lett. A 165, 409 (1992); J. D. 
Ramshaw, Phys. Lett. A 175, 169; L. Borland, A. R. 
Plastino and C. Tsallis, J. Math. Phys. 39, 6490 (1998). 

[13] C. Tsallis, Physica A 221, 277 (1995). 

[14] B. H. Lavenda and J. Dunning-Davies, Found. Phys. 
Lett. 3, 435 (1990); B. H. Lavenda, Statistical physics: 
A probabilistic approach (Wiley, New York, 1991) (and 
see also the references cited there) 

[15] K.T. Waldeer and H.M. Urbassek, Physica A 176, 325 
(1991). 

[16] P. T. Landsberg and R. B. Mann, Class. Quantum Grav. 

10, 2373 (1993) 
[17] C. Tsallis, R.S. Mendes and A.R. Plastino, Physica A 

261, 534 (1998). 
[18] P. J. Jund, S. G. Kim, and C. Tsallis, Phys. Rev. B 52, 50 

(1995); B. A. Mello et al. Phys. Lett. A 244, 277 (1998) 



4 



